Complete Angular Distribution Measurements of 
Two-Body Deuteron Photodisintegration between 0.5 and 3 GeV 
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Nearly complete angular distributions of the two-body deuteron photodisintegration differential 
cross section have been measured using the CLAS detector and the tagged photon beam at JLab. 
The data cover photon energies between 0.5 and 3.0 GeV and center-of-mass proton scattering angles 
10°-160°. The data show a persistent forward-backward angle asymmetry over the explored energy 
range, and are well-described by the non-perturbative Quark Gluon String Model. 
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I. INTRODUCTION 



Quantum Chromodynamics (QCD) has been success- 
fully applied in describing the structure and production 
of hadrons at high energies where perturbation theory 
can be used. There one can derive QCD scaling laws 
for the cross sections and hadronic helicity conservation 
laws. However, nuclear reactions have been convention- 
ally described in terms of baryons and mesons rather than 
quarks and gluons. It is therefore interesting and impor- 
tant to know in which energy region the transition from 
hadronic picture to quark-gluon picture takes place. This 
is why major efforts in nuclear physics have been de- 
voted, both theoretically and experimentally, to looking 
for qualitatively new phenomena which arise from the 
underlying quark degrees-of-freedom, and which cannot 
be modeled using meson field theories. 

Deuteron photodisintegration at high energies is es- 
pecially suited for this study, because a relatively large 
amount of momentum is transferred to the nucleons for 
a relatively low incident photon energy 0, • One pos- 
sible signature for the transition from nucleon-meson to 
quark-gluon degrees of freedom is the scaling of reaction 
cross sections above some incident photon energy. In par- 
ticular, simple Constituent Counting Rules (CCR) 0, H 
predict an asymptotic s^ 1 1 dependence of da/dt of the 
process at all proton angles. Here s and t are the invari- 
ant Mandelstam variables for the total energy squared 
and the four-momentum transfer squared, respectively. 

Deuteron photodisintegration cross-sections above 
1.2 GeV are available for photon energies E 1 up to 
5 GeV at three center-of-mass proton angles, $S m ' = 



36°, 52°, 69°; up to 4 GeV at tf£ m - =90° HE HSU; 
and at eight angles with i?£ m - = 30° - 143°, for E 1 = 1.6, 
1.9, and 2.4 GeV [T(J]. The asymptotic scaling predicted 
by CCR is observed at i?£ m - = 69° and 90° already at 
E 1 = 1 GeV and at $ c p m - = 52° and 36° only from 3 and 
4 GeV, respectively. In contrast, polarization observ- 
ables measured at t?p m ' = 90° for photon energies up to 
2 GeV ^1 0| do not support hadronic helicity conser- 
vation predicted by pQCD. Thus, it seems that although 
the observation of the scaling in the cross section at a few 
proton angles indicates the onset of the quark-gluon de- 
grees of freedom, the appropriate underlying physics has 
a mixture of perturbative and non-perturbative QCD as- 
pects. 

In this context, several non-pQCD models attempt 
to account for the experimental results using different 
strategies. The Reduced Nuclear Amplitude model 
(RNA) ^3| incorporates some of the soft physics not de- 
scribed by pQCD by using experimentally determined 
nucleon form factors to describe the gluon exchanges 
within the nucleons. The RNA calculation is only avail- 
able at -d c p m ' = 90° and makes no predictions for the an- 
gular dependence of the cross section. The calculations 
are normalized to data at energies sufficiently large, as- 
suming that perturbative regime is reached. 

The Hard Quark Rescattering Mechanism model 
(HRM) 0113 assumes that the photon is absorbed by 
a quark in one nucleon, followed by a high momentum 
transfer with a quark of the other nucleon leading to 
the production of two nucleons with high relative mo- 
mentum. The nuclear scattering amplitude is expressed 
as a convolution of the large pn scattering amplitude, 
the hard photon-quark interaction vertex and the low- 
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momentum nuclear wave function. The authors use ex- 
perimental data for the pn cross section, but since data 
do not exist for the actual kinematic conditions needed, 
they must be extrapolated, and predictions for deuteron 
photodisintegration are given as a band corresponding to 
the uncertainties introduced by the extrapolations. The 
model provides a parameter- free prediction of da/dt at 
i9p' m ' = 90°, and introduces a phenomenological function 
f(t/s) that is close to unity at flp m - = 90°, and varies 
slowly with Another attempt 0] to describe the 

deuteron photodisintegration within the same theoreti- 
cal framework of HRM, using an exact calculation of the 
quark exchange amplitude, provides evidence that the as- 
sumption used in [lj, [13 are questionable. 

The Quark Gluon String model (QGS) [13,111113 de- 
scribes the reaction as proceeding through three-quark 
exchange, with an arbitrary number of gluon exchanges. 
The exchanged nuclcon is replaced by a nucleon Regge 
trajectory that represents the sum of a tower of ex- 
changed nucleon resonances. The best description of the 
data is obtained using a nonlinear Regge trajectory. The 
model takes all but two of its free parameters from other 
processes, and fixes the remaining two using the experi- 
mental data on the deuteron photodisintegration cross 
section at E 1 = 1.6 GeV and i?£- m - = 36° and 52°. It 
provides the angular distributions and polarization ob- 
servables for few-GeV beam energies. 

Despite appearances, hard deuteron photodisintegra- 
tion is an intractable problem in meson-baryon theo- 
ries, the Asymptotic Meson Exchange Current model 
(AMEC) is able to extrapolate the conventional N— n 
interaction mechanisms to higher energy using form fac- 
tors to describe the dNN interaction vertex, and fix an 
overall normalization factor by fitting the experimental 
data at 1 GeV. 

A better insight into the competing models can be ob- 
tained from more detailed angular distributions of dif- 
ferential cross sections over broader angular and energy 
ranges than those presently available, and for final states 
involving different polarizations of the final hadrons. 

We report here the first measurement of nearly- 
complete angular distributions (10° < i} c p m - < 160°) of 
the two-body deuteron photodisintegration cross section 
obtained with the CEBAF Large Acceptance Spectro- 
meter (CLAS) in Hall B at the Thomas Jefferson Na- 
tional Accelerator Facility (experiment E93-017) [2l| for 
photon energies between 0.5 and 3 GeV. The data offer 
the opportunity for a detailed study of the energy de- 
pendence of differential cross section of the reaction at 
fixed proton angles, aiming at determining the onset of 
asymptotic scaling |22|. 

In the following, we first give some details of the experi- 
ment (Sec. II) and its data analysis (Sec. III). Then, we 
present our results on the deuteron photodisintegration 
cross sections da/dfl and da/dt, and compare them to 
available theoretical models and existing data (Sec. IV). 
We conclude with a summary (Sec. V). 



II. EXPERIMENTAL SETUP 

The data described in this paper were collected at the 
Thomas Jefferson National Accelerator Facility (JLab) 
during a 32-day run in August and September 1999 us- 
ing the Hall B tagged photon beam [23| and the CLAS 
detector [24j. The bremsstrahlung photon beam was 
produced by a 10 — 13 nA continuous electron beam of 
energy E = 2.5 GeV (August) and 3.1 GeV (September) 
impinging on a gold foil of 10 radiation lengths. A tag- 
ging spectrometer, with an energy resolution of O.li^o %, 
was used to tag ~ 10 7 photons per second in the energy 
range (0.20 - 0.95)E . 

A cylindrical mylar cryogenic target 10 cm long and 
4 cm in diameter was filled with liquid deuterium at 
about 23.7 K. The final-state particles were detected in 
the CLAS spectrometer, which is built around six super- 
conducting coils producing a toroidal magnetic field sym- 
metric about the beam and oriented primarily in the azi- 
muthal direction. The coils naturally separate the detec- 
tor into six sectors, each functioning as an independent 
magnetic spectrometer. Each sector is instrumented with 
3 sets of multi-wire drift chambers for track reconstruc- 
tion and one layer of scintillator counters, covering the 
angular range from 8° to 143° , for time-of-flight measure- 
ments. The forward region (8° < 1? < 45°) contains gas- 
filled threshold Cherenkov counters and lead-scintillator 
sandwich-type electromagnetic calorimeters for particle 
identification. For two CLAS sectors the coverage of the 
electromagnetic calorimeters is extended up to polar an- 
gles of 70°. The trigger for the data acquisition was de- 
fined by the coincidence between a signal in the tagger 
(identified photon) and one charged hadron in CLAS. Un- 
der these conditions 1771 million events were collected. 



III. DATA ANALYSIS 
A. Data selection 

A data quality check was performed to select runs with 
stable beam and detector performance. First, several 
run-based parameters normalized to the incident pho- 
ton flux were required to be constant at the few percent 
level from run to run: a) the total number of charged 
particles, and b) the number of particles p, 7r + , 7r _ , K + , 
and K~ . Then, c) the number of triggers with at least 
one charged particle in the final state for each tagger- 
timing counter; d) the number of photodisintegration 
events per 100 MeV; and e) the number of photodisin- 
tegration events per CLAS sector were required to be 
stable within the statistical errors. 

After applying the above data quality criteria, about 
7% of the originally collected data had been discarded. 
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B. Event selection 

Photodisintegration events — > pn were identified as 
follows: 

• The software coincidence time- window between the 
tagger and CLAS was set to ±1 ns, since the ma- 
chine electron bunches are separated by 2.004 ns. 

• Only events with a single charged particle, the 
proton, in the final state were selected. Pro- 
tons were identified by determining momentum and 
path length using the drift chambers, and velocity 
from the time-of-flight counters. 

• The reconstructed vertex position of the proton 
along with beam line was used to remove events 
originating outside the target cell. 

• Cuts on the square of the missing mass 
M x = (P 7 + Pd — P p ) 2 were performed to select 
exclusive two-body deuteron photodisintegration 
events. Here P 7 , Pd, and P p are the four-momenta 
of the photon, deuteron, and proton, respectively. 
In this study Mx is the mass of the neutron. 
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FIG. 1: Distributions for the Mx peak values for ■yd — * pX 
events before (top) and after (bottom) the momentum cor- 
rections are applied. The width of the corrected distribution 
is smaller (13.2 MeV rms — > 3.2 MeV rms) and the centroid 
is closer to the neutron rest mass. 



C. Momentum correction 

The momentum of the charged particles measured with 
CLAS strongly relies on the correct knowledge of the 
magnetic field geometry and the positioning of the drift 
chambers. Due to the complexity of the detector, and 
particularly of the superconducting torus magnet assem- 
bly, it is crucial to make sure that the momentum deter- 
mined by the drift chamber tracking system is reliable. 
For this reason the position of the peak of the miss- 
ing mass distributions from jd — » pX events has been 
checked over the whole range of proton momenta and 
scattering angles. After correcting for the energy loss in 
the target, the value of the peak was slightly off with re- 
spect to the neutron rest mass depending on the proton 
scattering angle. 

A correction procedure has been applied to the data 
using an empirical function depending only on the mea- 
sured three-momentum of the proton. It was assumed 
that the proton track angles are correctly measured, since 
the CLAS angular resolution is much better than the mo- 
mentum resolution [24[. We have also checked that the 
contribution due to the photon energy uncertainty is neg- 
ligible, by using exclusive jd — > ppix~ events. The cor- 
rection function has been calculated for each kinematic 
bin by fitting the ratio of the expected momentum, as 
calculated from the photon energy and theproton scat- 
tering angle, to the measured momentum |25| . 

The correction procedure introduced a significant im- 
provement in the resulting width and position of the peak 
in the missing mass distributions. Figurc^shows the dis- 
tributions for the peak values of the Mx distributions for 



7<i — > pX events for all CLAS sectors and all runs, before 
and after the correction was applied. Clearly, after the 
correction the distribution of the peak values is sharper 
(13.2 MeV rms ^3.2 MeV rms), and its mean value is 
closer to the neutron rest mass. 



D. Efficiency 

The single-proton detection efficiency in CLAS cannot 
be extracted from the deuteron photodisintegration data 
itself over the whole kinematic region of emitted pro- 
tons. Exclusive events, where both the neutron and the 
proton are detected, are limited because neutrons could 
be detected in the calorimeters only over a small angular 
range tf LAB < 45° for four CLAS sectors and tf LAB < 70° 
for the other two sectors. Other reaction channels with 
additional particles are problematic because the photo- 
disintegration protons have the highest momentum for a 
given proton angle and photon energy. 

For these reasons, the single proton efficiency has been 
evaluated using a GEANT simulation (GSIM) of the 
CLAS detector. Photodisintegration events have been 
generated uniformly in proton momentum and angle in 
the laboratory system, and then have been analyzed fol- 
lowing the standard reconstruction chain |2fi| . 

The proton detection efficiency egsim has been cal- 
culated in each kinematic bin in the laboratory system 
as the ratio of reconstructed protons -/Vrec to generated 
ones -/Vqen: 
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Bin widths for proton momentum AP^ —0.1 GeV/c 
and polar scattering angle Ad^ AB = 10° have been cho- 
sen. A smaller azimuthal angle bin of width A<^ LAB = 5° 
has been selected to better investigate the azimuthal be- 
havior of the CLAS proton detection efficiency, which 
gets worse on the boundaries of each sector due to the 
presence of the magnet coils. 

As an example, Pig.Elshows the resulting proton detec- 
tion efficiency for P p LAB = 0.95 GeV and t? lab = 65° as 
a function of the azimuthal angle <^ LAB . Similar plots 
have been obtained for the other proton angles. For 
^p AB = 45° — 125°, the proton efficiency is nearly con- 
stant in the central region of each sector, with an aver- 
age value of (94 ± 1)%, and drops sharply near the sector 
boundaries. At forward angles, the average efficiency de- 
creases dropping to about (50 ± 1)% at tf LAB = 15°. 




egsim are available with good statistics (proton momenta 
in the range 0.5 — 1.1 GeV and central regions of the six 
CLAS sectors) is shown in Fig. [21 The mean value of 
the distribution is very close to unity (0.997 ± 0.003). 

We checked also that regions of lower efficiency in 
CLAS corresponding to dead time-of-fiight paddles or 
drift chambers wires are well-reproduced by the simu- 
lation. Thus, the comparison over limited kinematics 
validates the GSIM results. 
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FIG. 3: Distribution of the ratio R — edata/egsim between 
the proton detection efficiency measured from the data using 
the jd — > pir~p reaction and that obtained from GSIM, for 
proton momenta in the interval 0.5 — 1.1 GeV and for the 
central 20° in ip for all sectors. 



FIG. 2: The CLAS proton detection efficiency evaluated 
using GSIM as a function of the azimuthal angle ip for 
P P LAB =0.95 GeV, i?p AB = 65°, and ip bins of 5°. 

In order to check the reliability of the simulations, the 
proton detection efficiency has also been obtained us- 
ing the data (where they are available) from the overde- 
termined jd — > ppn~ reaction. Each time a p and tt~ 
pair is found, the missing mass is calculated and the 
three-momentum is computed for candidates within tight 
constraints on the proton missing mass. The ratio be- 
tween the number of exclusive pir~p events, N p7r p , found 
by the particle identification, and the number of events 
with X identified as a proton by the missing mass cuts, 
N p7T x (p\ minus the number of the background events 
under the missing mass peak of the pTt~ X distribution, 
Nb, gives the experimental detection efficiency: 

£DATA = NP *-x(p) _ Ng ■ ( 2 ) 

The distribution of the values of the ratio 
R ■ = edata/cgsim calculated where both edata and 



E. Fiducial cuts and mean efficiencies 

As shown in Fig. [21 the proton detection efficiency is 
constant in the central azimuthal regions of the six CLAS 
sectors and decreases steeply near the sector boundaries. 
Thus, only events in a fiducial region (i.e. azimuthal re- 
gion of the phase space where the efficiency is uniform) of 
the detector have been used. For each bin in proton mo- 
mentum and scattering angle, and for each CLAS sector 
S, a mean efficiency e is defined as: 

e(AP p LAB ,A7? LAB ,5') = (3) 
<6 GSIM ( AP BAB , Atf LAB , S) ) ■ 77 ( AP p LAB , Atf LAB , S) . 

in which (e G siM (AP p LAB , A?? LAB , S)) is the average pro- 
ton detection efficiency over the fiducial Aip region, and 
ry (AP p LAB , Atf LAB , S) is the portion of the CLAS sector 
inside the fiducial cuts (i.e. the fraction of the (p interval 
considered). 
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F. Background subtraction 

At all proton angles and photon energies, the missing- 
mass distributions of the 7c? — > pX reaction show a neu- 
tron mass peak M n riding on a smooth background. As 
an example, figure0]shows the missing-mass distribution 
obtained for incident photon energy E~ = 0.95 GeV and 
proton scattering angle $ LAB = 25°. The logarithmic 
scale emphasizes the background contribution. These 
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result obtained from the extrapolation. The values have 
been found to be in a very good agreement with each 
other. 




M x (Ge\ 



FIG. 5: The behavior of the background contribution k as 
a function of the photon energy for proton scattering angles 
30° < tip m - < 40° and CLAS Sector 6. 



FIG. 4: Typical missing-mass spectrum of the reac- 
tion 'yd —> pX obtained in Sector 5 for photon energy of 
E-, =0.95 GeV and proton scattering angle i?p AB = 25°. 

distributions are well-reproduced with a Gaussian plus 
exponential form. Events within ±3er of the neutron peak 
have been kept for the determination of the cross section 
(here a is the width of the Gaussian distribution). 

The background contribution Nb to the number of 
total events under the peak N pea k has been evaluated 
by integrating the exponential fit function between the 
missing-mass cuts. At photon energies higher than 
(2.0 — 2.4) GeV depending on the proton angle, the neu- 
tron mass peaks are less clearly identifiable due to the 
low statistics. In these cases, the missing mass cuts for 
the selection of exclusive events have been obtained us- 
ing a second-order polynomial fit in $p m - and E 1 , of the 
(M P eak + 3<7) and (Mpoak — 3cr ) values determined at lower 
photon energies. In each $p m - bin the background has 
been evaluated by using a linear extrapolation of the fits 
in Ery to the ratio k — NB/N pca k obtained at lower pho- 
ton energies. As an example, Fig. [3] shows the values 
of the ratios k obtained for 30° < 0£ m - < 40° for CLAS 
Sector 6. Similar plots are obtained for other proton an- 
gles and CLAS sectors. Here k increases with the photon 
energy and proton angle $°- m - due to the loss in momen- 
tum resolution. 

In order to check the extrapolation procedure in the 
photon energy region above (2.0 — 2.4) GeV, the back- 
ground contribution has been evaluated from the data 
using larger bins (then increasing the statistics and mak- 
ing clearly identifiable the peaks) and compared to the 



G. Photon Flux 

The incident flux of photons on the target is given by 
N 7 = N e -ex, where N e is the number of tagged electrons, 
and ct is the tagging efficiency. N e has been measured 
online during the production runs, while ex has been 
measured during the normalization runs at low inten- 
sity (~ 10 5 7/s) using a nearly 100% efficiency lead-glass 
total absorption counter. We assume that ex remains un- 
changed during the production runs. Normalization runs 
were performed every time the experimental conditions 
for production runs were changed: a total of 178 million 
normalization events were collected. Fig. [B] shows the 
tagging efficiency measured in all the normalization runs 
at Eq~2.5 GeV. Similar results were obtained in the nor- 
malization runs at Eq~3.1 GeV. The tagging efficiency 
is stable at a level of ~ 2%. 



H. Systematic Uncertainties 

The contributions to the overall systematic uncertainty 
come from i) the determination of the number of incident 
photons, 1.9%, evaluated by looking at the variation 
of the number of photons per tagger channel in norma- 
lization runs; ii) the determination of the target length 
and density, w 0.5%; Hi) the proton detection efficiency, 
(2 - 8)%, evaluated as (eqsim - £data)/ £gsim; and iv) 
the background subtraction, around (1 — 2)% for photon 
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FIG. 6: Tagging efficiency for the 61 tagger-timing counters, 
measured during the Eq = 2.5 GeV normalization runs. In 
some runs (open points) the low energy counters had been 
switched off to increase the statistics at high photon energies. 



energies below 1 GeV and higher (up to ~ 6%) at for- 
ward and backward angles where the detector resolution 
and efficiency are worse. The latter has been evaluated 
by repeating the data analysis using both missing-mass 
cuts reduced and enlarged by 20%, and looking at the 
variation of the differential cross section. 

The resulting total systematic error is < 10% in the 
whole measured range. 



IV. RESULTS 

The photodisintegration cross section was calculated 
using: 



7' u p 



(4) 



pxN A N-y(E-y)An 



where N^, ak is the number of jd — > pn events weighted 
by the efficiency, Aft is the solid angle, A is the target 
molecular weight, Na is Avogadro's number, p is the tar- 
get density and x the target effective length. 

In Figures [7\ and [H] the angular distributions dcr/dVt 
are shown as a function of #p' m ' for photon energy bins 
100 MeV wide, in the range from 0.5 up to 3.0 GeV, 
and proton scattering angle bins 10° wide in the range 
10° < 0£ m - < 160°. The data are also given in Tables 
and [n] They are averaged over the six CLAS sectors. 
The results obtained by the six CLAS sectors separately 
are consistent with each other within the systematic er- 
rors. 



This is the first measurement of the nearly complete 
angular distributions of the — > pn reaction for pho- 
ton energies between 0.5 and 3.0 GeV. It allows one to 
investigate the behavior of the cross section in the very 
forward and backward angular regions. The data show 
a clear forward/backward angle asymmetry in the whole 
range of explored photon energies. At high energies the 
cross sections increase at very forward and backward an- 
gles. 

Also shown in Figures [7| and |S] are the previous data 
and the predictions of the few available models. For 
E 1 = 0.5 - 0.6 GeV, the Mainz data [H are slightly 
higher than CLAS at intermediate scattering angles. 
Starting from E 1 — 0.7 — 0.8 GeV the comparison can 
be extended also to the SLAC |E IS an d JLab Hall 
C data. The CLAS results agree well with the 
data from these experiments. For E 1 = 1.6 — 1.7 GeV, 
E y = 1.9 - 2.0 GeV, and E 1 = 2.4 - 2.5 GeV, the CLAS 
results agree with the angular distributions measured 
by the JLab Hall A collaboration 0] (the latter 
cover a smaller range in the proton scattering angle: 
26° < i?p' m ' < 143°), and extend to the very forward and 
backward angular regions where the cross section in- 
creases. 

For E~j > 1.0 GeV, the predictions of the QGS mo- 
del are shown in Figures [7| and |S] (solid curve). 
This model describes the angular distributions very well, 
and accounts for the persistent forward/backward angle 
asymmetry seen in the data by invoking the interference 
of the isovector and isoscalar components of the photon. 
The interference is constructive at forward angles and de- 
structive in the backward direction. The most forward 
points support the presence of the local maximum at 
about 20° predicted by the model. The backward points 
do not extend far enough to check for the presence of 
the second maximum. Also shown in Figures and |H1 are 
the predictions of the HRM model 27] (hatched band) 
calculated using the best angular fit for fixed energy pn 
scattering data. The bands reflect the poor accuracy of 
the data. The model agrees reasonably well with data in 
the central angular region over the whole explored energy 
range, and is lower at forward and backward angles apart 
from E 1 = 1.8 — 2.5 GeV. This agreement suggests a 
further investigation as, in principle, the HRM model is 
applicable for energies greater than ~ 2 GeV. 

The rich amount of CLAS data has made a detailed 
study of the power law dependence s~ n of the differen- 
tial cross section da/dt possible, in order to determine 
the onset threshold for the appearance of the s -11 scal- 
ing law predicted by perturbative QCD. This study 
indicates a proton transverse momentum scaling thre- 
sholds of Pt = 1.0 — 1.3 GeV/c for angles between 60° 
and 130°, and 0.6 — 0.9 GeV/c for forward and backward 
angles, with a nearly symmetric behavior around 90°. 

Figure [5] shows the results of da/dt (full circles) multi- 
plied by the factor s 11 predicted by CCR and plotted as 
a function of E y for the four proton scattering angles for 
which the predictions from all existing models are avail- 
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able. Also shown in the figure are the previous data: 
Mainz H (open squares), SLAC H S Q (full/green 
down-triangles), JLab Hall A 0] (full/blue squares) and 
Hall C 0,13 (full/black up-triangles). The two points at 
the same energy value from [lOj shown in the top panel 
come from two slightly different proton angles (30.3° and 
37.4°). The HRM model H3 (hatched band) agrees rea- 
sonably well with data up to about 4 GeV, then tends 
to be higher at forward angles. The RNA calculation 
is only available at #£- m - = 90°. The estimate for this 
figure |23| (dashed lines) is normalized to the datum at 
E 1 = 3.16 GeV. Other estimates at different angles have 
been provided in other papers by different authors but 
suffer from an incorrect normalization [29j, and there- 
fore are not shown in the figure. The AMEC model [2£j 
(dotted lines) predicts a slightly different energy depen- 
dence. The data at forward angles suggest a slower de- 
crease of the cross section with energy than predicted. 
Moreover there is a discrepancy for the highest energy at 
60° — 70° . Surprisingly, the model strongly overestimate 
data at energies lower than 1.6 — 2.0 GeV. The QGS mo- 
del describes well the data at all four proton angles. The 
largest discrepancy is found at 30° — 40° above 3 GeV 
where it suggests a slower decrease of the cross section 
with energy than observed. 

Clearly, further theoretical developments in this non- 
perturbative regime would be desirable to understand the 
transition region between the meson exchange picture 
and the QCD description of high energy nuclear reac- 
tions. 



Recherche Scientifique and the Commissariat a l'Energie 
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tional Science Foundation, and the Korea Science and 
Engineering Foundation. The Southeastern Universi- 
ties Research Association (SURA) operates the Thomas 
Jefferson National Accelerator Facility for the United 
States Department of Energy under contract DE-AC05- 
84ER40150. 



V. SUMMARY 



Differential cross sections for 7c? — > pn have been mea- 
sured for the first time with a nearly complete angular 
coverage (10° < i9p m - < 160°) in the photon energy range 
from 0.5 to 3.0 GeV using the CLAS detector and the 
tagged photon beam of Hall B at Jefferson Lab. The 
shapes of the angular distributions da/dfl show a persis- 
tent forward/backward angle asymmetry over the whole 
explored energy range. The cross sections da/dt fall by 
2-3 orders of magnitude from 1 to 3 GeV photon energy. 
The data have been used to determine the scaling thre- 
shold at every proton angle between 10° and 150°. The 
non-perturbative Hard Quark Rescattering Mechanism 
and Quark Gluon String models describe the data well. 
The latter accounts well also for the forward and back- 
ward angle asymmetry. 
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FIG. 7: (Color) Angular distributions of the deuteron photodisintegration cross section measured by the CLAS (full/red circles) 
in the incident photon energy range 0.50 — 1.70 GeV. Results from Mainz [2(jf (open squares, average of the measured values 
in the given photon energy intervals), SLAC 0-0 (full/green down-triangles), JLab Hall A [l0| (full/blue squares) and 
Hall C @, @ (full/black up-triangles) are also shown. Error bars represent the statistical uncertainties only. The solid line and 
the hatched area represent the predictions of the QGS fTsf and the HRM |27| models, respectively. 
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FIG. 8: (Color) Same as Fig.Qfor photon energies 1.7 - 3.0 GeV. 
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FIG. 9: (Color) Deuteron photodisintegration cross sections s dcr/dt as a function of Ej for the given proton scattering an gles . 
Results from CLAS (full/red circles), Mainz |2ql (open squares), SLAC 00 LI (full/green down-triangles), JLab Hall A 
(full/blue squares) and Hall C @,|9] (full/black up-triangles) are included, as well as predictions of the QGS (solid line), 
AMEC [20| and RNA [l3| models (dotted and dashed lines, respectively), and the HRM |2?| model (hatched area). Error bars 
represent the statistical uncertainties only. 
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± 0.049 


0.121 


± 0.121 


± 0.006 


0.072 


± 


0.072 


± 


0.003 












2.651 


0.257 


± 


0.182 


± 


0.016 


0.315 


± 


0.183 


± 


0.015 


0.183 


± 


0.134 


± 


0.009 


0.063 


± 


0.063 


± 


0.002 


2.749 


0.428 


± 


0.248 


± 


0.028 


0.237 


± 


0.168 


± 


0.011 












0.148 


± 


0.105 


± 


0.006 


2.886 


0.224 


± 0.225 


± 0.015 


0.134 


± 0.134 ± 0.006 


0.115 


± 


0.115 


± 


0.005 


0.110 


± 


0.110 


± 0.004 



TABLE I: Differential cross sections da/d^l in nb/sterad of the deuteron photodisintegration for photon energies 0.5-3.0 GeV 
and for center-of-mass proton angles 10° < m ' < 90° . The first error is statistical and the second is systematic in each case. 
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(E~<) (GeV) 


90° < 


e c p m - < ioo° 


100° < 


/lc.m. 
V p 


< 110° 


110° < 


/)C.m. 
O p 


< 120° 


120° < e c p m - 


< 130 


0.560 


595.8 


± 


5.1 


± 


23.8 


545.4 


± 


5.0 


± 


22.3 


477.9 


± 


6.9 


± 


25.9 


407.8 


± 


7.1 


± 


23.9 


0.641 


394.9 


± 


3.6 


± 


15.9 


358.6 


± 


3.1 


± 


14.7 


322.1 


± 


3.4 


± 


17.5 


258.0 


± 


3.5 


± 


15.2 


0.750 


209.3 


± 


2.8 


± 


8.6 


185.3 


± 


2.4 


± 


7.7 


171.2 


± 


2.6 


± 


9.4 


145.0 


± 


2.8 


± 


8.6 


0.864 


86.76 


± 


1.18 


± 


3.62 


85.01 


± 


1.10 


± 


3.57 


77.02 


± 


1.16 


± 


4.26 


57.99 


± 


1.18 


± 


3.49 


0.950 


50.64 


± 


0.72 


± 


2.16 


53.61 


± 


0.70 


± 


2.27 


47.33 


± 


0.73 


± 


2.64 


37.02 


± 


0.68 


± 


2.25 


1.052 


31.12 


± 


0.59 


± 


1.36 


30.12 


± 


0.59 


± 


1.29 


29.15 


± 


0.59 


± 


1.64 


21.77 


± 


0.58 


± 


1.35 


1.142 


17.30 


± 


0.41 


± 


0.77 


18.01 


± 


0.43 


± 


0.78 


17.09 


± 


0.42 


± 


0.97 


13.27 


± 


0.43 


± 


0.83 


1.254 


10.43 


± 


0.33 


± 


0.48 


10.55 


± 


0.34 


± 


0.47 


8.62 


± 


0.31 


± 


0.50 


7.49 


± 


0.33 


± 


0.48 


1.355 


6.36 


± 


0.26 


± 


0.30 


6.48 


± 


0.27 


± 


0.29 


5.23 


± 


0.24 


± 


0.31 


5.21 


± 


0.27 


± 


0.34 


1.449 


3.69 


± 


0.20 


± 


0.18 


3.79 


± 


0.20 


± 


0.17 


3.45 


± 


0.19 


± 


0.21 


3.17 


± 


0.21 


± 


0.21 


1.548 


3.13 


± 


0.19 


± 


0.16 


3.02 


± 


0.19 


± 


0.14 


2.17 


± 


0.16 


± 


0.13 


2.53 


± 


0.20 


± 


0.17 


1.648 


2.17 


± 


0.16 


± 


0.11 


1.83 


± 


0.15 


± 


0.09 


1.94 


± 


0.15 


± 


0.12 


1.97 


± 


0.17 


± 


0.14 


1.751 


1.36 


± 


0.13 


± 


0.07 


1.39 


± 


0.13 


± 


0.07 


1.27 


± 


0.13 


± 


0.08 


1.18 


± 


0.13 


± 


0.09 


1.845 


1.03 


± 


0.12 


± 


0.06 


0.694 


± 


0.099 


± 


0.034 


1.22 


± 


0.13 


± 


0.08 


1.11 


± 


0.14 


± 


0.08 


1.959 


0.633 


± 


0.109 


± 


0.037 


0.748 


± 


0.109 


± 


0.037 


0.599 


± 


0.103 


± 


0.039 


0.719 


± 


0.111 


± 


0.055 


2.051 


0.600 


± 


0.105 


± 


0.036 


0.643 


± 


0.100 


± 


0.032 


0.638 


± 


0.111 


± 


0.042 


0.456 


± 


0.086 


± 


0.036 


2.147 


0.310 


± 


0.109 


± 


0.019 


0.310 


± 


0.072 


± 


0.016 


0.479 


± 


0.134 


± 


0.032 


0.562 


± 


0.102 


± 


0.045 


2.250 


0.766 


± 


0.203 


± 


0.049 


0.304 


± 


0.076 


± 


0.016 


0.420 


± 


0.086 


± 


0.029 


0.284 


± 


0.072 


± 


0.024 


2.331 


0.297 


± 


0.094 


± 


0.019 


0.367 


± 


0.123 


± 


0.019 


0.243 


± 


0.115 


± 


0.017 


0.403 


± 


0.107 


± 


0.034 


2.458 


0.250 


± 


0.125 


± 


0.017 


0.116 


± 


0.083 


± 


0.006 


0.236 


± 


0.118 


± 


0.017 


0.117 


± 


0.083 


± 


0.010 


2.548 


0.134 


± 


0.095 


± 


0.009 


0.188 


± 


0.110 


± 


0.010 


0.058 


± 


0.058 


± 


0.004 


0.191 


± 


0.112 


± 


0.017 


2.651 


0.181 


± 


0.105 


± 


0.013 


0.056 


± 


0.056 


± 


0.003 


0.162 


± 


0.094 


± 


0.012 


0.214 


± 


0.108 


± 


0.020 


2.749 


0.139 


± 


0.098 


± 


0.010 












0.199 


± 


0.115 


± 


0.015 


0.123 


± 


0.087 


± 


0.012 


9 RR(\ 

i.OOU 


0.230 


± 


0.164 


± 


0.018 


0.066 


± 


n nfifi 


± 


0.004 


0.066 


± 


0.066 


± 


0.005 


0.063 


± 


0.063 


± 


0.006 


(E,) (GeV) 


130° < 


nc.m. 

P 


< 140° 


140° < 


oc.m. 
O p 


< 150° 


150° < e c p m - 


< 160 




0.560 


383.2 


± 


5.5 


± 


18.0 


347.4 


± 


12.1 


± 


15.4 


149.2 


± 


15.6 


± 


10.8 












0.641 


231.3 


± 


3.2 


± 


10.9 


207.7 


± 


4.0 


± 


9.2 


73.75 


± 


5.31 


± 


5.33 












0.750 


126.4 


± 


2.5 


± 


6.0 


120.7 


± 


3.2 


± 


5.4 


53.32 


± 


4.45 


± 


3.85 












0.864 


59.20 


± 


1.14 


± 


2.80 


56.80 


± 


1.37 


± 


2.58 


41.52 


± 


2.42 


± 


3.00 












0.950 


38.31 


± 


0.74 


± 


1.81 


39.69 


± 


0.89 


± 


1.81 


29.99 


± 


1.44 


± 


2.17 












1.052 


23.56 


± 


0.65 


± 


1.12 


27.65 


± 


0.82 


± 


1.27 


19.28 


± 


1.13 


± 


1.39 












1.142 


13.28 


± 


0.46 


± 


0.63 


18.36 


± 


0.61 


± 


0.85 


17.35 


± 


0.96 


± 


1.25 












1.254 


7.44 


± 


0.35 


± 


0.35 


14.58 


± 


0.55 


± 


0.68 


13.96 


± 


0.84 


± 


1.01 












1.355 


5.90 


± 


0.31 


± 


0.28 


12.09 


± 


0.49 


± 


0.57 


13.15 


± 


0.80 


± 


0.95 












1.449 


4.08 


± 


0.25 


± 


0.20 


9.34 


± 


0.42 


± 


0.45 


13.71 


± 


0.79 


± 


0.99 












1.548 


2.81 


± 


0.22 


± 


0.13 


7.29 


± 


0.39 


± 


0.35 


10.50 


± 


0.70 


± 


0.76 












1.648 


1.59 


± 


0.15 


± 


0.08 


4.66 


± 


0.30 


± 


0.23 


8.05 


± 


0.59 


± 


0.58 












1.751 


1.20 


± 


0.13 


± 


0.06 


3.33 


± 


0.26 


± 


0.16 


5.99 


± 


0.51 


± 


0.43 












1.845 


1.07 


± 


0.14 


± 


0.05 


2.55 


± 


0.24 


± 


0.13 


4.90 


± 


0.48 


± 


0.35 












1.959 


0.692 


± 


0.114 


± 


0.034 


1.15 


± 


0.16 


± 


0.06 


2.49 


± 


0.35 


± 


0.18 












2.051 


0.507 


± 


0.110 


± 


0.025 


1.33 


± 


0.17 


± 


0.07 


3.10 


± 


0.38 


± 


0.22 












2.147 


0.320 


± 


0.080 


± 


0.016 


0.942 


± 


0.146 


± 


0.048 


2.38 


± 


0.34 


± 


0.17 












2.250 


0.555 


± 


0.134 


± 


0.027 


0.552 


± 


0.113 


± 


0.028 


1.62 


± 


0.28 


± 


0.12 












2.331 


0.481 


± 


0.139 


± 


0.023 


0.531 


± 


0.161 


± 


0.027 


1.05 


± 


0.29 


± 


0.08 












2.458 


0.166 


± 


0.096 


± 


0.008 


0.364 


± 0.164 


± 


0.019 


0.832 


± 


0.343 


± 


0.06 












2.548 


0.327 


± 


0.135 


± 


0.016 


0.278 


± 0.140 


± 0.015 


1.08 


± 


0.39 


± 


0.08 












2.651 












0.310 


± 


0.140 


± 


0.017 


0.325 


± 


0.189 


± 


0.023 












2.749 


0.166 


± 


0.096 


± 


0.008 


0.136 


± 


0.097 


± 


0.007 


0.535 


± 


0.276 


± 0.039 












2.886 


0.064 


± 


0.047 


± 


0.003 


0.072 


± 0.073 


± 0.004 


0.516 


± 0.197 ± 0.037 













TABLE II: Differential cross sections da/dfl in nb/sterad of the deuteron photodisintegration for photon energies 0.5-3.0 GeV 
and for center-of-mass proton angles 90° < #p' m ' < 160°. The first error is statistical and the second is systematic in each case. 



